A note on bigravity and dark matter 
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We show that a class of bi-gravity theories contain solutions describing dark matter. A par- 
ticular member of this class is also shown to be equivalent to the Eddington-Born-Infeld gravity, 
recently proposed as a candidate for dark matter. Bigravity theories also have cosmological de Sitter 
backgrounds and we find solutions interpolating between matter and acceleration eras. 

PACS numbers: 



Observations show that most of the energy density in 
the universe is in the form of dark matter and dark en- 
ergy [U, H, H, 0, Hi ■ It is therefore of importance to have 
a simple and natural candidate for these components. 

Motivated by Yang-Mills theory, multi-graviton ac- 
tions are attractive extensions of general relativity. Con- 
sider several metrics g^^, a = l..n. The physical proper- 
ties and geometric interpretation of such a theory present 
great challenges. It is known that the full (diff)" sym- 
metry cannot be preserved by consistent interactions Q . 
The most general action preserving the full symmetry 
group is a sum of decoupled Einstein-Hilbert terms for 
each metric. 

However, interesting theories can be built by break- 
ing the (diff)" symmetry down to the diagonal subgroup. 
For n = 2, a particular "bi-gravity" theory with metrics 
{df^^^qt^A and action 



(1) 
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has been extensively studied 0, i, i, [M [O, E El, Q 
[l5| . Here A and A are cosmological constants for each sec- 
tor of the theory, g'"' and g^'^ are the inverses of q^i, and 
g^u respectively. K is the Ricci scalar for the metric g^jy. 
K is a dimensionless coupling. A class of bi-measure the- 
ories have been considered in (iGj and references therein. 

The interaction term above was first proposed in Ref. 
Q (see note [1^). A key feature of this interaction term 
is that it gives rise to Fierz-Pauli mass terms [17|] for the 
spin-2 fields. The above is not the most general mixing 
term that satisfies this condition. In particular, the den- 
sity -y/^ could be replaced by (—?)"(— 5) for any 
real u. However, for this theory to give rise to a dark 
matter dominated era, we find that, under the assump- 
tions described below, u — 1/2 is required. 

In this short note we point out the following proper- 
ties of ([1]). First we prove that for k = the action ([T|) 



is equivalent to the Eddington-Born-Infeld (FBI) theory 
proposed in as a theory for dark matter and dark 
energy. When one generalizes to the case k 7^ 0, it is 
a natural question to ask whether or not the dark mat- 
ter/dark energy interpretation still holds. The answer 
is in the affirmative. The metric q^i, can behave both 
as matter or as dark energy, and there exist solutions 
interpolating between them. We present two types of 
de Sitter vacua and study their stability under a certain 
set of perturbations. One is the well known solution in 
which the metrics are proportional. In the other case, 
which has received less attention, the de Sitter line el- 
ements of the two metrics are not proportional. These 
type of backgrounds were pointed out in different con- 
texts (e.g., in [1^ for the flat case, and in [lO] in static 
coordinates). We also display what are the conditions on 
the couplings that determine whether the Universe evolv- 
ing from a matter era ends up in the proportional or not 
proportional vacuum. Finally, we analyze tensor fluctua- 
tions on the de Sitter backgrounds. For the proportional 
case the equations can be decoupled and a condition on 
the couplings ensuring absence of tachyons is displayed. 

The large scale structure of the action ([1]) with k = 0, 
called FBI theory, has recently been studied in Ref. [2l[ . 
In that article, it was shown that the FBI theory has 
a phase for which the Friedmann background evolu- 
tion, growth of inhomogeneities and Cosmic Microwave 
Background (CMB) angular power spectrum are indistin- 
guishable from those predicted by ACDM. These results 
provide extra support for these theories as candidates for 
dark matter/dark energy. We shall come back to this 
point at the end. 

Theories interpolating between dark matter and dark 
energy are not new. Examples are the Chaplygin 
gas{2l m ill and the roUing tachyonjH, [H, 27]. For 
the Chaplygin gas, observational consistency of this in- 
terpolation has been challenged in Ref. [2^. 

We start by analyzing the relationship between the ac- 
tion ([T]) when K — and the FBI theory written in p^ . 
This is straightforward. In what follows it is convenient 
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to refer all lengths to I. We thus define two new dimen- 
sionless couplings a and ap by 



A 



On 



A 



a 



(2) 



The first step is to write the action ([T]) in Palatini 
form with a connection C^„^ associated to the metric 
Qfiv Varying with respect to g^i^, the equations can be 
algebraically solved for this field, 



K„ 



1 



(3) 



and therefore the solution can be replaced back in the 
action. The resulting action depends on g^j^ and the 
connection C^^a and is precisely the EBI action, 
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{R - 2A) + 



(4) 



Note that the action ([T]), for k = 0, in its Palatini form, is 
a parent action in the sense that one may either eliminate 
Qfiu to obtain EBI, or eliminate the connections to get ([T]) 
in terms of the metrics only. 

Born-Infcld type actions have appeared repeatedly for 
many years in many different contexts. The action (jj]) 
is particularly close to the one discussed in [2§] although 
different in interpretation. Another class of class of Born- 
Infeld theories are tachyonic fields [1^, [2^, which are an- 
other candidate for dark matter and dark Energv[30|. 
This field is described by the following effective action, 



It=- 



du. 



4>) + lEHig^ip), 

(5) 

where lEH{9fj.u) is the Einstein-Hilbert action including a 
cosmological constant and V((j)) is an effective potential, 
which, in open string theory is|3l| 



Vo 



cosh(a(/)) 



(6) 



Just as the square root in the EBI theory ([4]) can be 
transmuted into a standard kinetic term by introducing 
g^y, a similar manipulation holds for ([5|). Consider the 
following action, 

■(7) 

This action can be seen as Polyakov's version of (O , ex- 
cept that the metric field is q^^. This is seen by varying 
with respect to q"^ . We obtain an equation which allows 
us to algebraically solve qap in terms of 0, therefore, we 
may put it back in the action ([7]). We get 
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(8) 



which is precisely the tachyon action ([5]) action when 
U{^) = -Vo/V{cl)). 

Considering the form of action fT]), it is suggestive to 
add a kinetic term to the auxiliary metric field q^^. The 
obvious choice is to also add a Einstein Hilbert term with 
a cosmological constant. If we do so, we obtain precisely 
the first line of ([T]) plus Ip in This is EBI action in 
bi-gravity form plus a scalar field minimally coupled to 
the metric q^^, . The variation of this action with respect 
to q^jy gives, again, an equation that may be solved al- 
gebraically for the g-metric. Inserting this back in the 
action, and redefining Vo, we obtain the following gener- 
alization of (HI, 



I 
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(9) 



^- det(g^i. - PKf,^ + df,(t)du(, 



which becomes EBI when U — (f) — 0. 

Our second goal is to study the cosmological properties 
of the bi-gravity system described by ^ . Assuming that 
both metrics are homogeneous and isotropic flat FRW, 



ds^ = -dt^ 



ds^ ^ -X^dt^ 



+ Y^dx^. (10) 



For a given set of couplings a,aQ, k there exists more than 
one de Sitter vacua. We distinguish two cases: propor- 
tional vacuums (PV) if both metrics are proportional and 
non-proportional vacuums (NPV). In the proportional 
case the functions a, X,Y are given by 



^2 ^ 1 - (qq - 6k) 



with 



1 — (ao — 6k)q; 
3(1 -a) ■ 



Y^aX (11) 



(12) 



The only condition for the existence of this solution is the 
positivity of the constants and X^ above. Note that 
if ao — 6k = ct~^, then this vacuum becomes Minkowski. 
Also note that for a close enough to 1, or k sufficiently 
large, the de Sitter acceleration can be made arbitrarily 
large, even if the cosmological constant, ao/P, vanishes. 
Conversely, even for big values of ao we may fine-tune 
the couplings in order to obtain arbitrarily small acceler- 
ation. This is an attractive feature in the context of the 
problem of the cosmological constant. 

A second class of de Sitter vacua with non-proportional 
metrics also exists (NPV). Let Y — aXA, where A is a 
constant (which is 1 for the previous case). In this case 
a — eT"* where A, X and H are constants determined by 
the equations: 



X^A^ 



ao = 

1 



X^A^ 



a 



4X2^2 ^2 



4X2 



(13) 
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To find the metric parameters H, A, X one needs to solve 
a third order algebraic equation. This means that, in gen- 
eral, we may expect three different NPV for a given set 
of couplings. In [T9^, these kind of solutions are also dis- 
cussed. In that case, however, the cosmological constants 
are adjusted to have flat backgrounds, so the mixing term 
used here, called Vi in that reference, gives rise only to 
the proportional vacua. 

We may ask now if the above vacua are stable. Con- 
sider perturbations of the form 

a{t) = a„{t) + ea^{t), 
X{t) - X„{t) + ^X,{t), 
Y(t) = Y^{t)+eY,{t), 

where the subscript ^ indicates the background solutions 
found above. We expand the equations to linear order in 
e and look for the conditions on the couplings such that 
the perturbations do not grow in time. These conditions 
are best expressed with a picture. In Fig. 1, we show 
the regions of stability for the PV and the NPV with 
ao = 1/2 and a — 0.9 respectively. 

The vertical axis represent the value of a. Note that 
there are regions where both solutions are stable. Those 
are the regions where we have degenerated vacuum. We 
do not know, however, if one of them turns out to be 
metastable. 
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I. The Friedmann equation for a(€) in the limit of a(t) 
sufhciently small reads. 



3ij2 K, SttGp 



^0 



1 



PXq a3 ■ 



(15) 



where p is the energy density of other conventional fluids 
(e.g. radiation or baryons). Therefore, we see that the 
q-metric plays the role of an additional dust-like matter, 
irrespective of the presence of other conventional fluids 
like radiation or baryons. Adjusting the constant 
one can have any desired amount of "dark matter" . In 

particular choosing — 3.34 x 10~'wcMpc~ where 
Wc ~ 0.09 — 0.12 we get the right amount of dark matter 
as required by cosmological observations. Note that the 
K constant has no role at the above regime, and so both 
the bi-gravity action ^ and the EBI model describes the 
same physics when a{t) <C 1. 

Our final task is to study gravitons propagating on the 
de Sitter vacua. In the following we use conformal time. 
We start by perturbing the FRW metrics as 



and 



-dr"^ + (7ij -|- hij)dx^dx-''\ 



(16) 



Y\J,,+X^J)dx'dx^ (17) 



where hij is the tensor mode perturbation of the (^-metric 
and Xij the tensor mode perturbation of the g-metric. 
The tensor modes are transverse and traceless. 

From now on we drop the indices on hij and Xij since 
no confusion arises. We find that the field equations for 
h and x are 

h + 2-h^^^h = 
a 

2a2 ,,^2 



[X^ - 2(1 - wo)^] [h - x) (18) 



FIG. 1: Stability regions for the proportional vacuum solution 
(PV) and the non-proportional ones (NPV). In the left the 
value of ao = 0.5 is fixed. In the right a = 0.9 is fixed. 



Now we proceed to show, by studying the background 
cosmological evolution, that this theory encompass a 
dark matter candidate, even for n ^ Q. To this end, the 
equations of motion coming from ([1]) are evaluated with 
homogeneous and isotropic metrics. We assume that 
both metrics are not singular at the same time (37j , which 
implies that, when the scale factor a{t) [with a(0) = 0] is 
sufficiently close to zero, 

X{t) « + Xxt, Y{t) « Vo + Yit, (14) 

where Xq^Yq, Xi,Yi are constants. Inserting this into 
the equations of motion, with q and g metrics given by 
(ITUl) . one finds Yi and Xi as a function of Xq, Iq, a, and 



where wo — 



and 




2a^uio 
'~PX^ 



wqV X 



- 2(1 - ii;o)«] {h-x) 



(19) 



respectively. We now adopt the above equations to the 
special case of de Sitter vaccua. For the proportional de 
Sitter vacuum described above with 4- = 0, = - and 
Wo = — 1 these two equations can be collected in matrix 
from as 



d d 



a dr 

where the mass matrix A^^ is 



-1 
1 



(20) 



(21) 



4 



One of the eigenvalues is clearly zero while the other one 
is 



(22) 



Thus in order for the theory not to contain spin-2 
tachyons, we must have X"^ > An which translates to 



1 - ao 



2(1 + 2a)K 



1 



> 



(23) 



If these conditions are fulfilled, this theory describes a 
massless and massive graviton. In particular, they imply 
stability for the EBI vacuum tensorial modes, which were 
first studied in Ref-lSa]. 

It is important to mention that gravitons are not the 
only physical excitations; vector and scalars modes may 
also propagate. The reason is that the action has two 
metrics but only the diagonal subgroup of diffeomor- 
phisms leaves the action invariant. This means that the 
scalar and vector modes of only one of the metrics can 
be set to zero by a gauge symmetry. This raises the issue 
of the stability of the theory which should be analyzed 
along the lines of Ref . [ssl . [34] . It was shown in those ref- 
erences that cosmological massive gravitons are stable if 
the mass of the graviton obeys m? > 2 A/3. We expect 
similar results to hold in our case. 

To conclude, we have studied in this note several cos- 
mological aspects of bigravity actions of the form ([1]). 



Most importantly we have shown that generically these 
actions contain a phase at early times where the second 
metric behaves as dark matter. The equations also ad- 
mit a de-Sitter background which is an attractor if some 
conditions on the couplings are fulfilled. This implies a 
transition between the "matter" and "de Sitter" phases. 
For the case k = this transition has been explored in 
detailed in [2lj and shown to be problematic at the level 
of fluctuations and the calculation of CMB spectra. How- 
ever, as shown in [2lj, one can choose initial conditions 
and couplings such that the metric q^^ is locked into its 
matter phase up until today. In this case, bigravity pre- 
dicts a CMB spectrum which is indistinguishable from 
standard particulate dark matter. The calculations for 
K ^ are far more complicated and we shall consider 
them in a separate publication. 
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